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m ' ABSTRACT 



We study the propagation of radio waves through a medium possessing density 
fluctuations that are elongated along the ambient magnetic field and described by an 



anisotropic Goldreich-Sridhar power spectrum. We derive general formulas for the 
wave phase structure function D§, visibility, angular broadening, diffraction-pattern 
length scales, and scintillation time scale for arbitrary distributions of turbulence 
along the line of sight, and specialize these formulas to idealized cases. In general, 
D§ (8r) 5 / 3 when the baseline 8r is in the inertial range of the turbulent density 
spectrum, and D§ °c (8r) 2 when Sr is in the dissipation range, just as for an isotropic 
Kolmogorov spectrum of fluctuations. When the density structures that dominate the 
scattering have an axial ratio R ^> 1 (typically R ~ 10 ), the axial ratio of the broad- 
ened image of a point source in the standard Markov approximation is at most ~ R l l 2 , 
and this maximum value is obtained in the unrealistic case that the scattering medium 
is confined to a thin screen in which the magnetic field has a single direction. If the 
projection of the magnetic field within the screen onto the plane of the sky rotates 
through an angle A\|/ along the line of sight from one side of the screen to the other, 
and if R ~ 1//2 <C A\\i <C 1 , then the axial ratio of the resulting broadened image of a point 
source is 2(8/3) 3 / 5 /Ay ~ 3.6/Ay. The error in this formula increases with A\|/, but 
reaches only ~ 15% when A\|/ = n. This indicates that a moderate amount of variation 
in the direction of the magnetic field along the line of sight can dramatically decrease 
the anisotropy of a broadened image. When R ^> 1, the observed anisotropy will in 
general be determined by the degree of variation of the field direction along the sight 
line and not by the degree of density anisotropy. Although this makes it difficult to 
determine observationally the degree of anisotropy in interstellar density fluctuations, 
observed anisotropics in broadened images provide general support for anisotropic 
models of interstellar turbulence. Regions in which the angle y between the magnetic 
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field and line of sight is small cause enhanced scattering due to the increased coher- 
ence of density structures along the line of sight. In the exceedingly rare and probably 
unrealized case that scattering is dominated by regions in which y < (Sr//) 1 / 3 , where 
/ is the outer scale (stirring scale) of the turbulence, D§ (8r) 4 / 3 for 8r in the inertial 
range. In a companion paper (Backer & Chandran) we discuss the semi-annual modu- 
lation in the scintillation time of a nearby pulsar for which the field-direction variation 
along the line of sight is expected to be moderately small. 



1. Introduction 

Scattering of radio waves from point sources by electron-density fluctuations in the interstel- 
lar medium (ISM) gives rise to a number of effects, including intensity scintillation and angular 
broadening (Rickett 1990). These phenomena provide a useful diagnostic of density fluctuations 
in the ISM on scales ~ 10 8 — 10 10 cm (diffractive scintillation) and 10 13 — 10 15 cm (refractive scin- 
tillation) (Armstrong et al. 1995). In a majority of cases, broadened images are anisotropic, with 
axial ratios (long dimension of image divided by short dimension) between 1.1 and 1.8 for sight 
lines through strongly scattering regions in the Galactic disk (Mutel & Lestrade 1990, Wilkinson 
et al. 1994, Molnar et al. 1995, Spangler & Cordes 1998, Trotter et al. 1998) and as large as 3:1 
for OH masers near the Galactic center (van Langevelde et al. 1992, Frail et al. 1994). [Axial 
ratios > 5 : 1 have been observed in the solar wind (Narayan et al. 1990, Armstrong et al. 1990).] 
This paper explores the relation between anisotropic scattering and recent theories of anisotropic 
magnetohydrodynamic (MHD) turbulence. 

Early theories of MHD turbulence assumed isotropy. Iroshnikov (1963) and Kraichnan (1965) 
independently derived a fc~ 3 / 2 power spectrum for the velocity and magnetic field in isotropic in- 
compressible turbulence. Over the last decade, however, a number of authors have investigated the- 
ories in which small-scale fluctuations are elongated along the local direction of the magnetic field, 
B (e.g., Montgomery & Turner 1981, Shebalin et al. 1983, Higdon 1984, Higdon 1986, Oughton 
et al. 1994, Sridhar & Goldreich 1994, Goldreich & Sridhar 1995, Montgomery & Matthaeus 
1995, Ghosh & Goldstein 1997, Goldreich & Sridhar 1997, Matthaeus et al. 1998, Spangler 1999, 
Bhattacharjee & Ng 2000, Cho & Vishniac 2000, Maron & Goldreich 2001, Lithwick & Goldreich 
2001). Using phenomenological arguments and a statistical turbulence theory, Goldreich & Sridhar 
(1995) (hereafter GS) derived a form for the velocity and magnetic power spectrum in incompress- 
ible MHD turbulence that approximately corresponds to the anisotropic spectra found in direct 
numerical simulations (Cho & Vishniac 2000, Maron & Goldreich 2001). Lithwick & Goldreich 
(2001) extended the GS theory to compressible MHD turbulence. They found that the spectra of 
shear-Alfven modes, slow modes, and entropy modes all have the same anisotropic form as the 
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shear-Alfven and pseudo-Alfven modes in incompressible turbulence. The density fluctuations in 
the compressible turbulence are dominated by the slow modes and entropy modes. Lithwick & 
Goldreich (2001) considered the effects of damping by neutrals, radiative cooling, electron heat 
conduction, and ion diffusion, and found that the density spectrum can extend to the small scales 
responsible for diffractive scintillation provided that the neutral fraction is very small, and that 
either |3 (the ratio of thermal to magnetic pressure) is not much larger than 1 or the outer scale / of 
the turbulence is fairly small. For |3 ~ 1, the density spectrum is cut off at wavelengths along the 
magnetic field comparable to the proton mean free path. Because the wavelength of the density 
fluctuations across the magnetic field is much smaller than the parallel wavelength, the density 
spectrum extends to an inner scale d significantly smaller than the proton mean free path, of order 
(Lithwick & Goldreich 2001) 

The axial ratio (long dimension divided by short dimension) of fluctuations of perpendicular scale 
d in the GS spectrum is ~ (Z/d) 1 / 3 , which is ~ 10 3 for / = 1 pc, (3 = 1, and n = 1 cm~ 3 . 

In this paper, we take the interstellar density fluctuations to have a GS spectrum and calculate 
the consequences for radio wave propagation in the ISM. In section 2 we review the relations be- 
tween the wave phase structure function, visibility, angular broadening, diffraction-pattern length 
scales, and scintillation time scale. In section 3 we derive formulas for these quantities for arbi- 
trary distributions of turbulence along the line of sight, and specialize these formulas to idealized 
observational scenarios. In section 4 we give a detailed summary of our main results. 



2. Background: wave phase structure function, visibility, angular broadening, 
diffraction-pattern length scales, and scintillation time scale 

In this section, some general results on scintillation and angular broadening are reviewed. For 
an overview of the subject, the reader is referred to Rickett (1990). A systematic derivation of the 
interferometric visibility and intensity correlation function for a plane wave propagating through 
a stationary medium given a stationary observer was given by Lee & Jokipii (1975a,b). Their 
derivation is extended to the case of a moving point source and moving observer in appendix A. 

The visibility, (E (r, t)E* (r + 8r, t + 8t) ) , is the correlation between the electric field observed 
at position r in the earth's reference frame at time t and the electric field at position f + 6> in 
the earth's reference at time t + 8t. In the Markov approximation (Lee & Jokipii 1975a), which 
requires that as a wave propagates through one correlation length of an electron density fluctuation 
the change to the wave field induced by the density fluctuation is small, the visibility is given by 
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(Lotova & Chashei 1981, Cordes & Rickett 1998; see appendix A for a derivation), 

(E(r,t)E*(r+8r,t + bt)) = cxp[-D^(br,bt)/2}, (2) 

where 

D^(br,6t) =4%r 2 e l 2 f dzg(z,Z(t)r,5t)) (3) 
J 

is the wave phase structure function, 

dq x / dq y { 

-OO — oo 

[1 - cos(q • a)] P„ e 9y, ft = 0; z) } , (4) 
d(?,8t) = j$r + V eg 8t, (5) 

and 

Veff(z)= (l-^)v p + ^V obs . (6) 

Here, P ne (q;z) is the power spectrum of the electron density fluctuations, q is the Fourier-space 
wave vector, z is the coordinate along the path from the source to the observer, L is the distance 
between source and observer, r e is the classical radius of the electron, X is the wavelength at which 
the observations are taken, V p and V bs are the velocities of the pulsar and observer, respectively, 
relative to the frame of the density fluctuations (which we assume to have negligible motion) and 

— * 

perpendicular to the line of sight, and V e ff is the effective perpendicular velocity of the sight line 
relative to the plasma turbulence. The power spectrum P„ e is taken to depend upon z in what is 
essentially a two-scale approximation, since the outer scale of the turbulence, /, is taken to be 
much smaller than L. The sight line sweeps across a small-scale density fluctuation in a time 
short compared to the Lagrangian correlation time of the density fluctuation, and thus the intrinsic 
time evolution of the density fluctuations is ignored. The wave phase structure function D§ can 
be thought of as the mean-square difference between the density-fluctuation-induced phase incre- 
ments along one line of sight from the source's position at time t — (L/c) to the observing location 
at time t and another line of sight from the source's position at t + 8t — (L/ c) to the observing loca- 
tion at time t + dt. For strong scattering [equivalent to D§ 3> 1 as Sr — > oo], the intensity correlation 

— * 

function < I(r,t)I(r+8r,t + 8t) > is approximately given by (Lotova & Chashei 1981, Cordes & 
Rickett 1998; see appendix A for a derivation) 

<I(r,t)I(r+$r t t + tit) >~< / > 2 {1 +exp[-D ()) (5r,8?)}, (7) 
where < / > is the mean intensity. 

The observed electric field and intensity vary stochastically as the sight line moves through 
turbulent density fluctuations in the ISM. To obtain the visibility and intensity correlation function, 
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the observed electric field and intensity are averaged over some interval in time. In the derivation 
of the formulas quoted above, it is assumed for simplicity that this averaging procedure is equiv- 
alent to taking an ensemble average over the electron density fluctuations which determine the 
observed electromagnetic fields. This assumption is justified only if the integration time is suffi- 
ciently long (Goodman & Narayan 1989). The effects of scattering by anisotropic turbulence for 
shorter-duration observations are beyond the scope of this paper. 

We define orthogonal coordinates x and y in the plane perpendicular to the line of sight such 
that in the case of anisotropic scattering, y is the direction of strongest scattering (largest angular 
size). The length scales of the diffraction pattern along the x and y directions, given by l± x and /(ir- 
respectively, and the scintillation time scale t d are determined from D§(Sx,8y, ht) = D§(hr,ht) by 
the equations 

Af(/«M,0,0) = 1, (8) 
D^(0,/ dj ,0) = 1, and (9) 
D^(0,0,t d ) = 1. (10) 

The effective angular size of the source in the x and y directions is given by (Rickett 1990) 

e.s,x = rr-, and (11) 



1 



®s,y = —■ (12) 



3. Calculation of the wave phase structure function, diffraction-pattern length scales, 
scintillation time scale, and angular broadening for a GS spectrum of density 

fluctuations 

We define orthogonal x, y, and z axes with z along the line of sight and y along the direction of 
maximum angular broadening. We define orthogonal x 1 , y', and z' axes with z' along z and with an 
angle between x' and x and between y' and y, as in figure 1 . The value of \\f will vary along the 
line of sight so that x' is aligned with the projection of B(z) in the plane of the sky. We also define 
orthogonal x" , y", and z" such that y" is along / in the original xy plane, and z" is parallel to B(z), 
which is taken to make an angle y(z) with the z and z' axes, as depicted in figure 2. The density 
fluctuations are elongated along B(z) and thus z". The density structures and scattering have the 
same statistical properties for J = Jo and y = n — Jq. The separation d between the two lines of 
sight in the phase- structure-function formula is taken to make an angle £(z) with respect to the xf 
axis, as in figure 1 . 
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Fig. 1. — The eccentric ellipse along the x' axis corresponds to the diffraction pattern associated 
with a thin layer of material with a single magnetic field vector whose projection in the plane of the 
sky, B^, is along x' . The projection of the density structures would also appear as eccentric ellipses 
elongated along x' . The less eccentric ellipse corresponds to the diffraction pattern resulting from 
scattering along the entire line of sight. The y axis is defined to be along the short axis of the less 
eccentric ellipse, which is the long dimension of the corresponding broadened image. 



We assume a GS spectrum for P ne (q,z) with a sharp cutoff at scales larger than an outer scale 
/ and at scales smaller than an inner scale d, i.e., if < q± < d~ l then 

P ne (q,z) = (l/6n)(dnl)r^q- l0/3 f {j^Iy^j > (") 



and if q± < I 1 or qj_ > d 1 then P Me (q, z) = 0, where 



I 1 if \x\ <C 1 

/Ws, !o itw>i' (14) 



and where q^ = (q^„ + qy,) 1 / 2 . The normalization in equation (13) has been chosen so that the 
mean square density fluctuation is given by (8ng). 
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Fig. 2. — The magnetic field is along z" , which makes an angle y with the z and z' axes. 
Noting that 

Qx" = <?x' cosy- <? z > sin y, (15) 
qyi = qy, and (16) 
Qz" = <? x 'siny+c7 z 'CosY, (17) 

(18) 

we find from equations (13) and (14) that P ne (q x r,qy,q z / = 0,z) is nonzero only when 

l^y/sinyl 3 < q 2 x d 2 cos 2 y+qy I 2 . (19) 

The fluctuations that dominate radio-wave scattering typically satisfy qyl ^> 1. When <jy/ ^> 1 and 
y^> (qyl)~ 1 / 3 (i.e., the magnetic field is not too closely aligned with the line of sight), the first 
term on the right-hand side of equation (19) is negligible compared to the second for all allowable 
values of q x /, and the upper limit on q x > becomes 



^'/siny| 3 <q 2 yl 2 . 



(20) 
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Thus, when qyl 3> 1 and y 3> (qyl) 1 / 3 , 

Pn e (qx>,qyi,q z ' = o,z) ^ 

(1 ^ )<5 "« >rl/ ^ 10/3/ (^fc Y )' 

On the other hand, as y ^ 0, P Me becomes isotropic in and <gy. 

In order to calculate D§ for different cases, we will first calculate the value of g(z,o) for 
various orientations of the magnetic field direction z" relative to the line of sight and the vector a 
separating the two lines of sight in the turbulent medium. 



3.1. Value of g(z, 6) when the magnetic field is along the line of sight 

When o is in the inertial range (d <^ o <^ I) and y <C (a//) 1 / 3 , the dominant contributions 
to g(z,d) comes from values of q x r and qy satisfying q 2 + q 2 , ~ a 2 , for which the value of / 
in P ne (q x ',q y ',<2V = 0,z) is equal to 1. Thus, when a is in the inertial range (d <C o <C /) and 
y< (o//) 1 / 3 , one can take 

p ne (q^qy^ = o,z) ^ (i/67i)(5n e 2 )r 1 / 3 (^+^)- 5/3 , (22) 

which implies that to lowest order in (a//) 

« fc3) = 2^§3)« rl/V/3 ' <23) 

When a is in the dissipation range (a <C J) and y <C (J//) 1 / 3 , one finds that to lowest order in 
(a/J) and (d/l) 

^(z,6) = (l/8)(5n2)/- 1 /3 J -2/3 a 2_ (24) 

The values of g in equations (23) and (24) are large compared to those that arise from other ori- 
entations of B due to the increased coherence length of the elongated density structures along the 
line of sight. 



3.2. Value of g(z, 6) when magnetic field is not along the line of sight and o is not along x ; 

Let 

u = qyl, (25) 
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GCOSC 

a= — p5 < 1, 

and let u m i n be some value of u satisfying 

l««min« \b\~ l . 



(26) 
(27) 

(28) 



— 1/3 —1/3 

For y^> u ■' and for u > u m [ n , P„ e is given by equation (21). Thus, for y^> u ■ ' , the contribution 



to g(z, o) from u > u m [ n is 

g(z,a)= e/ / rfww 10/3 

•'"mill 



w 2//3 cscy- 



- sin(<2M 2//3 esc y) cos(bu) 

If J < osinC, < / and |<a| 3 / 2 (csc y) 3 / 2 <C which implies that 

I^IXa/O^^cscy) 3 / 2 , 



(29) 



(30) 



2/3 > 

then there exists a u\ < l/d such that au^ cscy^C 1 and \bui\ ^> 1. The contribution to g(z,o) 
from <7y 3> /~ ! is then dominated by values of u in the interval (w m i n , «i), in which the integrand 
in equation (29) can be expanded: 



g(z,o) = 2 ^" e ^ I duu 8 / 3 cscy[l -cos(^w)]. 

371 Ju.„-,„ 



(3D 



To lowest order the limits of integration can be changed to (0, °°), which implies that the contribu- 
tion to g(z,o) from values of qy much greater than is 



5r(5/3) 



(6n 2 )o 5 / 3 /- 2 / 3 cscy| sin^| 5 / 3 

if d < osin^ < I. 

mm 



(32) 



Although equation (32) was derived assuming y 3> w~ in , it is approximately valid under the 
slightly more general condition y 3> (c| sin^l/7) 1 / 3 since g(z,a) is dominated by values of 
of order (o\ sin^|) _1 . 

As C, — > 0, the contribution to d) from large qy vanishes, but the contribution from small 
qy does not. Anticipating the results of the next subsection, we note that in order for the large-gy 
contribution to dominate and for equation (32) to be accurate when osin^ is in the inertial range, 
C, must satisfy the condition 

KI>(o//) 1/5 (siny) 3 / 5 (33) 
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as well as equation (30). 

When a sin £ is in the dissipation range («C d) and y^$> (d/l) 1 / 3 , one finds that the contribution 
to g(z, a) from <jy 3> is to lowest order 

g(z,S) = (l/7r)(8^)/- 2 /3j-l/3 a 2 sin 2^ cscy 

if a sin ^< J. (34) 

In order for the large-gy contribution to dominate over the small-gy contribution, C, must be 3> 

(J//) 1 /6( s iny) 1 /2. 



3.3. Value of g(z, 6) when magnetic field is not along the line of sight but d is along x' 

When osin^ is in the inertial range, y > (csin^/7) 1 / 3 , and \Q -C (o//) 1 / 5 (siny) 3 / 5 , the dom- 
inant contribution to g(z, d) in the Markov approximation comes from qy ~ giving in order of 
magnitude 

g(z,o)~(&i2)o 2 r 1 . (35) 

Equation (35) also holds in the Markov approximation when o sin C, is in the dissipation range, y ^> 
(J//) 1 / 3 , and \Q <C (J//) 1 / 6 (siny) 1//2 . The Markov approximation, however, is formally invalid 
for density fluctuations at scales ~ /, where I is on the order of 1 to 100 pc for typical interstellar 
conditions. This is because the Markov approximation assumes that 

8(|) = r e Xl c 8n e <^ 1, (36) 

where 8(|) is the phase increment (relative to propagation through the average medium) induced 
as the wave propagates through one correlation length l c of the density fluctuation 8n e . For the 
Goldreich-Sridhar spectrum with fractional density fluctuations of order unity at the outer scale 
/, one has 8n e ~ n e (l c /l) 1 ^. Thus, the Markov approximation is valid only for those density 
fluctuations with correlation lengths satisfying 

/ c «/ 1 / 4 X- 3 / 4 r e - 3/ V 3/4 . (37) 

For / = 10 20 cm, X = 30 cm, and n e = 0.1 cm , the Markov approximation is only valid for 
l c 10 14 cm. Although the Markov approximation breaks down for the fluctuations at scales 
larger than ~ 10 14 cm, we will for simplicity assume that the Markov-approximation results are 
accurate for the large-scale fluctuations in order of magnitude. In almost all cases, this assumption 
leads to the conclusion that the large-scale density fluctuations are unimportant for scintillation 
and angular broadening. The only exception is when all of the density fluctuations are aligned 
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in approximately the same direction, as in section 3.6. In that case, the large-scale fluctuations 
dominate the angular broadening in the direction parallel to the elongated density structures. 

Having determined the values of g(z, d) for different orientations of the magnetic field relative 
to z and a, we now turn to calculations of D§. 



3.4. Value of when there is either significant variation in the magnetic-field direction z" 
along the line of sight, or when z" is not along the line-of-sight and x' is not along a 

The z integral appearing in equation (3) can be divided into intervals in which (a) the magnetic 
field B is along the line of sight, (b) B is not along the line of sight and x' is not along d, and (c) 
B is not along the line of sight but x' is along o. In intervals of type (a), when o is in the inertial 
range, the value of the integrand in equation (3) is ~ (//a) 1 / 3 times larger than in intervals of 
type (b). However, intervals of type (a) [for which y< (a//) 1 / 3 ] are about a factor of (a//) 2 / 3 
less common than intervals of type (b) when the field direction varies significantly along the line 
of sight in a random manner, and thus intervals of type (a) do not contribute significantly to D§. 
Intervals of type (c) are much less common than intervals of type (b), and in addition the integrand 
in equation (3) is smaller in intervals of type (c) than in intervals of type (b). Thus, D§ is dominated 
by intervals of type (b), both when there is significant random variation in the field direction along 
the line of sight and when the field direction is fixed in a single direction of type (b). 

When d <C a <C /, for most values of C, one also has d <^ osinC, <^ I, the condition under 
which g(z, d) is given by equation (32). Thus, when o is in the inertial range, 

A» = 5f^3)^ 2 / L Jz((5n ' )r2/3cSCYl sin ^ 5/3oV3 ) 

ifd<o<Z, (38) 

where the integration excludes segments of the line of sight for which y < (a//) 1 / 3 . Similar argu- 
ments show that when o is in the dissipation range, the z integral in equation (3) is dominated by 
intervals in which g(z, 3) is given by equation (34), giving 

D^ = Ar 2 e l 2 I ^((S^r^-^Vsin^cscy) 

J 

if a < d, (39) 
where the integration excludes segments of the line of sight for which y < (d/l) 1 ^. 

A subtlety not addressed in this paper is that o can be within the inertial range of the turbulence 
over part of the line of sight and within the dissipation range over the remainder of the line of 
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sight. For simplicity it will be assumed that if a is in the inertial range at either the source or at the 
observer, then the inertial-range formulas can be used for the entire line of sight. For example, if 
8r is in the inertial range and dt = 0, then o is taken to be in the inertial range along the entire line 
of sight, although in fact it is in the dissipation range sufficiently close to the source. 

From equations (8), (9), (10) and (38), one finds that for d <C a /, 

k,x= ^^^Lrig.cscYlsin^l 5 / 3 ^^ ' ^ (40) 

k,y= ^^^Lrigscscylcos^l 5 / 3 ^ 5 / 3 ^ and (41) 



td 



^L %s cscY|sinC| 5 / 3 |(l-^)y p + ^ obs | 5 / 3 | 



1 -3/5 

(42) 



where r| gs = [87t/5r(5/3)]r;;?i 2 (8ng)/~ 2//3 , £, = z/L, and the integration excludes segments of the 
line of sight in which y< (a//) 1 / 3 . Equations (40) and (41) show that the length scale of the 
diffraction pattern tends to be dominated by fluctuations near the observer. If V p 3> Vobs, then 
equation (42) shows that the time scale of the diffraction pattern tends to be dominated by fluctua- 
tions near the source. 



3.5. Special case: homogeneous extended medium with a statistically isotropic random 

magnetic field direction 

Let us suppose that the magnetic-field direction varies randomly in an isotropic manner along 
the line of sight, and that (Sn 2 .), I, and d are constant along the line of sight. We have already 
assumed that / <C L, which implies that in the z integral in equation (38) we can average over the 
direction of B while holding z constant to a good approximation. With 

327l 2 f n l 2 

a i 



[ ^(sinQ 5 / 3 , (43) 
Jo 



5r(5/3) Jo 
, and 

H = (r 2 l 2 (Snl)L)-\ (45) 



% = and (44) 



one finds that 



D§ = aH- l r 2 / 3 J o d\ | (6> - V p dt + VobsSO^ + V p 5r| 5 / 3 

if d < o < /. (46) 
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When dt = 0, or V p = V obs = 0, 

D^ = (3/S)aH- l r 2/3 8r 5 / 3 . if d <Sr<Z. (47) 
The diffraction pattern in this case is isotropic, with 

3/5 



/ 8 \ ' 

ha = k,y = ( 3^ ) HV5iV5 if d « V = V « /• 



(48) 



When o is in the dissipation range, one finds that 

D^ = 4n 2 H- l r 2 / 3 d- l/3 x 
[(I /3)\& + V ohs 8t-V p 8t\ 2 + br + V ohs 8t 2 } ifo<J, (49) 

(50) 



Ida = ky = (J^J ' H l ' 2 l l l 3 d l l 6 if h, x = l dj « d, 
3 \ 1/2 H l / 2 l l / 3 d 1 / 6 



and 



4 *V (Vl + Vobs-Vp + Vp 2 ) 1 / 2 

if/d,x = /d, y <^ (51) 



3.6. Special case: one thin scattering screen with a single direction of B 

In this section, it is assumed that all of the electron density fluctuations lie within a screen 
of thickness Az <C / <C L. The magnetic field direction thus does not rotate very much within the 
screen, and the density fluctuations are aligned in the same direction. The screen is taken to be a 
distance z from the source. Upon defining 

h = (r 2 e X 2 hz{hnl))- 1 (52) 

(for X = 10 cm, Az = 10 18 cm, and (Sn 2 ) = 1 cm" 6 , h = 1.3 x 10 5 cm), one finds that 

/d ' y= (^lr^) 3/5/z3/5/2/5(sinY)3/5 7 Xd^zld,y/L^l. (53) 
To obtain an approximate value for l A x , one must use equation (35) to obtain 

/d.^^ 1/2 / 1/2 -. (54) 

z 
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Equation (54) holds whether zld,x/L is in the dissipation range or the inertial range. As discussed 
in section 3.3, we are unable to determine how ld, x changes with y as esc y becomes large. When 
esc y is of order unity (e.g., < 4) the axial ratio of the anisotropic diffraction pattern and broadened 
image is given by 

!r* = ~ frY /10 Xd«zl d , y /L«l, (55) 
oj,x ld,y \nj 

ignoring factors of order unity. The axial ratio of the density structures in the Goldreich-Sridhar 
spectrum that dominate the diffraction pattern is 

f I \ 1/5 

[l/{zk,y/L)\ V 3 ~ I - J if J«z/ dj /L«/. (56) 

Thus, for cscy of order unity, the axial ratio of the diffraction pattern is roughly the square root 
of the axial ratio of the turbulent structures that dominate the scattering at the scale of the diffrac- 
tion pattern (i.e. turbulent structures with <jy 1 ~ zld,y/L)- It should be noted, however, that this 
result follows from applying the Markov approximation to the outer-scale density fluctuations, 
for which the approximation is formally invalid. There exist other anisotropic power spectra for 
which the axial ratio of the diffraction pattern arising from a thin screen equals the axial ratio of 
the density structures. For example, if P ne (lx"^ ( ly"^ ( lz") 06 + + {R") 2 q 2 ,]~^^ 2 where R" and 
|3 are constants, /?"» 1, and (3 < 4, and if cscy is of order 1, then the axial ratio of the diffrac- 
tion pattern is ~ R" (Backer & Chandran 2000). Similarly, Narayan & Hubbard (1988) showed 
that if the two-dimensional power spectrum of the phase fluctuations in a thin scattering screen is 
oc [{R") 2 q 2 + q 2 ]^^ 2 with P < 4, then the axial ratio of the associated diffraction pattern is R" . The 
key difference between these axial ratios and the axial ratios for a GS spectrum is the following: 
although the image broadening along the y direction for a GS spectrum is dominated by highly 
anisotropic small-scale fluctuations, the image broadening along the x direction for a GS spectrum 
is dominated by the large-scale fluctuations. Since the image broadening along x is enhanced rel- 
ative to the amount of image broadening arising solely from the small-scale density structures, the 
image is less anisotropic than the density structures. 

The scintillation time scale td is given by 



td 



5r(5/3) 1 3/5 /? 3 / 5 / 2 / 5 |cscC|(siny) 3 / 5 
8tt J 

if</«zVefffd/£«/- (57) 



However, as £ — > 0, td does not approach °o but instead reaches a maximum value when D§ is 
determined using equation (35): 

td,msx. ~ —r 7 • (58) 

Veff 
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If |V e ff is constant but C, varies continuously through a 360° interval during the course of either 
the earth's orbit around the sun or the orbit of a binary pulsar, then the ratio of the maximum to 
minimum scintillation times when esc y is of order unity is given by 



td.max ( I 



1/10 



(59) 



For zld,y/L in the dissipation range, one finds that 

h y = (l/2)/j 1 / 2 /V3ji/6( siny )i/2_ ifd / L< j. ( 60 ) 

z 

Since equation (54) applies in the dissipation range, the axial ratio of the diffraction pattern and 
broadened image for esc y of order unity is given by 

®s,y _ h,x ( ' 



9.s,x ld, : 



f-J if^ d , y /L«J, (61) 



which again, is approximately the square root of the axial ratio of the turbulent structures at scale 
d that dominate the scattering for zld, y /L <C d. The scintillation time scale is 

/*l/2/l/3jl/6( siny) l/2 csc £ 

^ = 2v7^ z dj/ <<c ( } 

As ^ — > 0, ^ stops increasing after reaching the value given in equation (58). Again, if V e ff is 
constant but C, varies continuously through a 360° interval during the course of either the earth's 
orbit around the sun or the orbit of a binary pulsar, then the ratio between the maximum and 
minimum scintillation times when esc y is of order unity is given by 



t, /n 1/6 



^d.min V d 



(63) 



3.7. Special case: thin scattering screen in which the direction of magnetic field in plane of 

sky rotates through an angle A\\f 

In this section, it is assumed that the scattering medium extends from z = z\ to z = z\ + Az, 
with Az L. It is assumed that (5«g) , /, and y axe constant within the screen, that o is in the inertial 
range of the density fluctuations, and that y^> (osin^//) 1 / 3 . The angle \\f between the projection 
of the magnetic field in the plane of the sky and the x axis is given by 

2 Az 
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The value of Ay is conceptually equivalent to Az/l. It is assumed that 



(65) 



where 7? is the axial ratio of the small-scale density structures that dominate the scattering. This 
means that for any direction of d, portions of the line of sight of type (c) in section 3.4 (for which 
d is nearly along the direction of B in the plane of the sky) contribute a negligible amount to D§, 
and D§ is given by equation (38). From this it follows that 



k,x 
k,y 



M\|//2 

/ d\|/|cos\|/| 5 / 3 
Jo 



A\|//2 



L JO 



d\\f\ sin\|/| 



5/3 



3/5 



(66) 



When \ A\\f \ <C 1, equation (66) reduces to 

l d , x ^ 2(8/3) 3 /5 3.6 



/do 



A\|/ A\|/' 



(67) 



Equations (66) and (67) are plotted in figure 3. If |V e ff| is constant during one orbit of the earth 
around the Sun or one orbit of a binary pulsar, then the ratio of the maximum to minimum scintil- 
lation times ?d.maxAd,min during the course of one orbit is also given by 



td,v 



Ay/2 



1 cos\|/| 5//3 



1 3/5 



A\|//2 



d\\f\ sinx^l 5 / 3 



(68) 



and when \A\\f \ < 1, 



^d,i 



2(8/3) 



3/5 



36 
A\|/' 



(69) 



^d,min A\)/ 

The scintillation time reaches its maximum value when the sight line moves parallel to x , the 
direction in which the density fluctuations are most elongated, and reaches its minimum when the 
sight line moves along y. 



4. Summary of results 



In this paper, we assume a GS spectrum (Goldreich & Sridhar 1995) of electron density fluc- 
tuations in the ISM. We then calculate general formulas for the wave phase structure function, 
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Fig. 3. — Axial ratio of the diffraction pattern and broadened image of a point source when the 
magnetic field direction in the plane of the sky rotates through an angle A\\f within the scattering 
screen. The solid line represents equation (66), and the dashed line represents equation (67), which 
is derived assuming \ A\\f \ <C 1. 

visibility, angular broadening, diffraction-pattern length scales, and scintillation time scale for ar- 
bitrary distributions of turbulence along the line of sight, and specialize these formulas to idealized 
cases. Our main results are as follows: 

(1) Unless the magnetic field is closely aligned with the line of sight over a much greater 
fraction of the line of sight than is expected for random magnetic fields, the scaling of the wave 
phase structure function with baseline 8r is D§ oc (8r) 5 / 3 for 8r in the inertial range of the 
turbulence, and D§ (Sr) 2 for Sr smaller than the dissipation scale or inner scale of the turbulence, 
just as for an isotropic Kolmogorov spectrum of density fluctuations. 

(2) If the magnetic field were closely aligned with the line of sight [within an angle ~ (5r/7) */ 3 
if Sr is in the inertial range, and ~ (J//) 1 / 3 if Sr is in the dissipation range] over more than a 
minimum-threshold fraction of the line of sight, scattering would be dominated by that fraction 
of the line of sight and D§ would be proportional to (Sr) 4 / 3 if 8r were in the inertial range, and 
to (8r) 2 if Sr were smaller than the dissipation scale. The minimum-threshold fraction is small 
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[~ (Sr//) 1 / 3 if 8r is in the inertial range, or ~ (J//) 1 / 3 if 8r is in the dissipation range], but much 
larger than expected for randomly varying magnetic fields. The increase in scattering when B 
becomes aligned along the line of sight is due to the increased coherence length of the elongated 
density structures along the line of sight. 

(3) Angular broadening is most anisotropic when scattering occurs within a thin screen in 
which the magnetic field has a single direction. In this case, the axial ratio of the broadened image 
of a point source in the Markov approximation is approximately the square root of the axial ratio 
of the density fluctuations that dominate the scattering, provided that cscy is of order unity (e.g. 
< 4), where y is the angle between the magnetic field and the line of sight. We have not found 
the dependence of the axial ratio on y when cscy is large, although we do show that if y = the 
broadened image is isotropic. 

(4) If the scattering medium is confined to a homogeneous thin screen of thickness Az <^ L, 
and if the direction of the magnetic field in the plane of the sky rotates through an angle A\|/ 3> 
R ~ l l 2 along the line of sight from one side of the screen to the other, where R is the axial ratio of 
the density fluctuations that dominate the scattering, then the axial ratio of the observed image of 
a point source is ~ 3.6/Ay. [This formula is derived assuming A\|/ <C 1, but closely approximates 
(to within 15%) the more general formula, equation (66), for A\\f as large as 7t.] This indicates 
that even a moderate variation in the field direction within the scattering medium can dramatically 
reduce the axial ratio of an angularly broadened imaged. If the field direction along the line of 
sight to a source varies significantly, the resulting angular-broadening anisotropy is determined by 
the amount of variation in the field direction rather than by the degree of anisotropy of the density 
fluctuations. 

(5) When the magnetic field in the scattering medium has a single direction or rotates through 
a relatively small angle A\|/, the scintillation time is longest when the sight line moves along the 
direction of greatest elongation of the electron density fluctuations in the scattering medium, and 
shortest when the sight line moves orthogonal to the direction of greatest elongation of the density 
structures. This gives rise to an observable modulation of the scintillation time during the course 
of either the Earth's orbit around the Sun, or of a binary pulsar's orbit. When the speed of the sight 
line in the plane of the sky V e ff in the scattering medium (assumed to be a thin screen) is constant 
during the orbit, the ratio of the maximum to minimum scintillation times *d,maxAd,min is the same 
as the axial ratio of the broadened image of a point source, when the magnetic field has a single 
direction in the thin screen as well as when the magnetic field rotates through an angle A\\r. If the 
field direction along the line of sight to a source varies significantly, the degree of variation in the 
scintillation time is determined by the amount of variation in the field direction rather than by the 
degree of anisotropy of the density fluctuations. 

(6) The diffraction-pattern length scales and angular broadening are more sensitive to fluctua- 
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tions near the earth than to fluctuations near the pulsar. The time scale for diffractive scintillations 
is more sensitive to fluctuations near the pulsar if the pulsar speed significantly exceeds the speed 
of the observer with respect to the interstellar turbulence. 

We thank Steve Spangler, Bob Mutel, Sridhar Seshadri, Jason Maron, and Barney Rickett 
for helpful discussions. This research has been supported by NSF grant AST-9820662 at UC 
Berkeley, and by NSF grant AST-0098086 and DOE grants DE-FG02-01ER54658 and DE-FC02- 
01ER54651 at the University of Iowa. 



A. Derivation of equations for visibility and intensity correlation function for a moving 

observer and moving point source 

In this appendix we derive equations (2) through (7) for the visibility and intensity correla- 
tion function for a moving point source and moving observer. In several places, the derivation is 
identical to the corresponding calculation for plane waves incident upon a turbulent medium and 
observer that are at rest, and the reader will be referred to Lee & Jokipii (1975a,b) for the details. 

We work in the rest frame of the source, and start with the scalar wave equation for electro- 
magnetic waves propagating in a cold non-magnetized plasma (Faraday rotation is ignored), 

, 1 d 2 E (Op 
^- 72W -^E = »- (AD 

where 

iol = — . (A2) 

' m e 

We look for solutions that are approximately spherical waves, 

E = ue i{ ~ kr - m \ (A3) 

with 

u — > 1 as r — > 0, (A4) 

where (r, 9, (j)) are spherical coordinates centered on the source, and, in the notation of Lee & 
Jokipii (1975a), 

9 G) 2 £o , . ro p 

fc=— ^, £ = <£*), and £* = 1 £. (A5) 

c z Or 

The angled brackets denote an ensemble average over the turbulent density fluctuations. We also 

define 

£ = -^ (A6) 
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to be the fractional fluctuations in the dielectric constant e*, where r e is the classical electron radius 

e 2 1 (m e c 2 ), so that 

£* =£ (l + e). 



We make the "parabolic approximations," 

du d 2 u 



, du () H 

k— > ^—r and a>— > 

dr dr 2 dt dt 2 



(A7) 



(A8) 



and using equations (Al) through (A8) find that 



du 1 di 



where 
and 



2ik\ ^ + — — + V1 w + e& z w = 0, 
dr v g dt 1 



v g = Cy/fy, 



-L 2 

r L 



i a 
sineae 



. n d\ i d 2 

sinG— + 



ae 



sin' 



ea^ 2 



We then transform coordinates from (r, 0, <|), t) to (z, 0, (|), s), where 

z = r and s = r — v g t. 



(A9) 

(A10) 
(All) 

(A12) 



Points with a constant 5 are all those values of r and t associated with the wavefront that arrives at 
the observer at the time t = (L — s)/v g , where L is the distance from source to observer. In terms 
of z and s and the shorthand notation defined in table 1 , equation (A9) can be written 



dU\ T T 

2ik— + ,ui+£ik z ui =0. 
dz 



(A13) 



Equation (A13) is of the same form as equation (7) of Lee & Jokipii (1975a). Taking the complex 
conjugate of equation (A13), and evaluating the terms at 02, <|>2, and S2, one can write 



du*2 

"a7 



2ik^ - V 2 - e 2 k 2 ul = 0. 



(A14) 



Multiplying equation (A 13) by u\ and equation (A 14) by u\ and adding, taking the ensemble aver- 
age, and noting that for statistically homogeneous density fluctuations V 2 ^wiwl) = ^\ 2(^1^2)' 
we find that 



2ik^-(u\u* 2 ) +k 2 ((£i -£ 2 )u\u* 2 ) = 0. 
dz 



(A15) 



We now make the Markov approximation, in which we assume that the correlation length of the 
density fluctuations / e is much smaller than the length /„ along the line of sight over which u 
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changes significantly. We then follow a procedure analogous to the one described in appendix A 
of Lee & Jokipii (1975a) to show that 

(( ei -e 2 ) MlM *) = ^|) f z £fe 1 ((eiel> + (e2^>-<ei4>-(e2ei». (A16) 

L J — oo 



To evaluate the right-hand side of equation (A 16), we first define a according to the equation 

(A17) 

We restrict our attention to values of a and § satisfying 
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6= 2 + - 



loci < 1 and 



<1. 



(A18) 



We then introduce Cartesian coordinates (x, y) in the plane perpendicular to the line of sight, with 



x = zOL and y = zty. 



(A19) 



(Note: the xyz coordinates used here are not related to and (|) in the conventional way. For 
example, the z axis corresponds not to = 0, but instead to = %/2 and (|) = 0.) We also make 
a two scale approximation along the line of sight, dividing the dependence of the variables on the 
radial coordinate into two parts: a dependence on Z reflecting variations over scales ~ l £ , and a 
dependence on z reflecting variations over scales / e . The integral in equation (A 16) effectively 
extends over a distance ~ l £ <C z, and so in calculating x and y for e[ and e' 2 we can use z in place 
of z' in equation (A 19) with only small error. Defining x = xx + yy + Zz, and assuming statistical 
homogeneity and stationarity, we can write 



($n e (xi,ti;z)5n e (x2,t2;z)) = Q(xi -x 2 ,h -h\z). 



(A20) 



Abbreviation 



Meaning 



V 2 



sin0 ; 30, 



u(z,Qi,tyi,Si) 
u(z',$u§usi) 
e(z, Bi,tyi,si) 
£(z',Bi,<S?i,Si) 

sin©^ ) + 



i a 2 

sin 2 ; a^) 2 



Table 1 : Shorthand notation, where i = 1 , 2, 3, or 4. 
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We assume that the density fluctuations are static in the rest frame of the turbulent medium, which 
is reasonable since the Lagrangian correlation time of the density fluctuations that dominate scin- 
tillation is small in the GS theory compared to the time for the line of sight to sweep across such a 
fluctuation. This gives 

Q(x,t;z) = G(x-U m t;z), (A21) 

where U m is the uniform velocity of the turbulent medium with respect to the source, and G(x;z) is 
the spatial correlation function of the density structures in the rest frame of the turbulent medium, 
which is the inverse Fourier transform of the density power spectrum P He : 



Defining 



we find that 



where 



G(x;z) = J d\S rx P n Xq\z). (A22) 
Hij= f dz'ieie'j + eje'i), (A23) 

J — oo 

Hij = 3 ~^- J d 3 qP ne (q;z)e r ^d (q 



Z V 8 



(A24) 



pij = z(a,- - a.j)x+z(to - tyj)y+ - Sj)Um . (A25) 



V 8 



The presence of the U m /v g term in the delta function in equation (A24) is due to aberration — in 
the frame of reference of the turbulent fluctuations, the wave appears to be moving in a direction 
slightly offset from the line of sight. Since U m v g , this effect is small and can be ignored, and 
the delta function can be written 8(k z ). Equation (A24) thus reduces to 



where 



J Jdq x dq y P ne (q x ,qy,q z = 0;z)e iUi \ (A26) 



Oij = z(ai - aj)x+z(to - $j)y+ - Sj)Um ^ 7 {All) 

Vn 



and where U m ± is the component of U m perpendicular to the line of sight 
Equations (A15), (A16), (A23), and (A26) imply that 



^{urnl) = -2nr 2 e l 2 g(z,Oi 2 )(uiul). (A28) 
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If the correlated observations are taken at locations that are separated in the observer's rest frame 

— * 

by a displacement 8r = 8xx + 8yy + 8zz and at times that are separated by a time interval 8t, then 

L(ai-a 2 ) = bx + U ohS7X 8t, (A29) 
1(01 -02) = 5y + t/ obs , y 5?, and (A30) 
Sl -s 2 = -v g (t 1 -t 2 ). (A31) 

(A32) 

Because u changes much more slowly along the line of sight than across the line of sight, 8z 
and £/ bs.z can be neglected. Defining V p and V b s as the velocities of the source and observer 
perpendicular to the line of sight measured in the rest frame of the turbulent medium, one has 

U m ,± = -V p , (A33) 

U ohs ,A- = -V p + V ohs , (A34) 

and 

012 = £(5r + ? ob8 Sf) + (l - |) V p dt. (A35) 
Integrating equation (A28) and using equations (A4) and (A35), one obtains equation (2). 

To derive the intensity correlation function, we define 

T2,2 = (U1U2U3I1I). (A36) 
Using the same procedure used to derive equation (A15), one obtains 

ar 

2ik^- + (Vi j + V 2 2 - V 2 3 - V 2 4 )r 2 , 2 + k 2 ((e l +£ 2 - e 3 - e 4 )uiu 2 u* 3 u$) = 0. (A37) 
oz ' ' 

Using the Markov approximation, we find that 

( (£l + £ 2 - £3 " e 4 )MlM2«3"4) = (2#1 1 + #12 + #34 - #13 - #14 - #23 - #24) , (A38) 

which gives 

d r 2, 2 i f„ 2 , vr2 vr2 x~72 "\ i -1 i ^2,2 

To an excellent approximation, 

V,.^ (,— + —). (A40) 



= ^ (Vi 5l + V 2 )2 - V 2 3 - V 2 )4 ) r 2 , 2 + — ^ (2ffi 1 + #12 + #34 - #13 " #14 " #23 ~ #2 4 ) • 

(A39) 



1 / a 2 a 2 

u ~ z 2 \da 2 a^ 2 



and thus equation (A39) is equivalent to equation (8) of Lee & Jokipii (1975b), if the effects of 
aberration are ignored. Lee & Jokipii's (1975b) derivation of the approximate intensity correlation 
function in the limit of strong scattering can therefore be applied to equation (A39) to obtain 
equation (7). 
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